Spin-orbit effects on two-electron states in nanowhisker double quantum dots 
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We investigate tlieoretically the combined effects of the electron-electron and the Rashba spin- 
orbit interactions on two electrons confined in quasi-one-dimensional AllnSb-based double quantum 
dots. We calculate the two-electron wave functions and explore the interplay between these two 
interactions on the energy levels and the spin of the states. The energy spectrum as a function 
of an applied magnetic field shows crossings and anticrossings between triplet and singlet states, 
associated with level mixing induced by the spin-orbit coupling. We find that the fields at which 
these crossings occur can be naturally controlled by the interdot barrier width, which controls the 
exchange integral in the structure. 
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I. INTRODUCTION 

Substantial efforts have been devoted to understand- 
ing and manipulating the electron spin aiming at poten- 
tial applications in semiconductor spintronics.^'^'"^ While 
the spin-orbit interaction in quantum dots has been 
extensively studied from the single-electron perspec- 
tive, its combined effects with the Coulomb interaction 
in few-particle systems has only recently begun to be 
exploredi^'^'^'i Especially in nanowhiskers, Fasth et al^ 
measured the strength of the spin-orbit interaction in 
two-electron InAs cylindrical dots (diameter ~ 50 nm, 
length ^ 120 nm), and Pfund et al^ studied spin re- 
laxation in a similar system. Here we examine two- 
electrons in a double quantum doiji^'iiii^'i^ system with 
only one transverse quantum mode active. The one- 
electron spectrum of this thin nanowhisker semiconduc- 
tor structure has been treated by us previously^^ as well 
as the phonon-mediated spin-relaxationii^ 

In this paper we calculate and analyze in detail the 
two-electron states in such a system, taking into account 
the Rashba spin-orbit coupling and the Coulomb inter- 
action between the electrons. We pay special attention 
to the degree of admixture of different two-electron spin 
wave-functions, which will influence the spin-flip transi- 
tions in this system. 

The paper is organized as follows. In Sec. [Til we mtro- 
duce the effective quasi-one-dimensional Hamiltonian of 
two interacting electrons in the presence of the Rashba 
interaction and describe the method and approximations 
used in the calculations. In Sec. Ull] we investigate the 
effects of the electron-electron interaction on the energy 
levels with and without Rashba interaction and their de- 
pendence on an applied magnetic field. We monitor the 
mean value of the spin projection as a function of the 
structural parameter that determines the strength of the 
Rashba spin-orbit coupling. In Sec. IIVI we provide some 
concluding remarks. 



II. THEORETICAL DESCRIPTION 

We investigate the problem of two interacting electrons 
in a quasi-one-dimensional double quantum dot structure 
in the presence of the structural or Rashba spin-orbit in- 
teraction. We study two identical 30 nm wide dots, sep- 
arated by an interdot barrier, 3 nm or 5 nm wide. In 
our calculations we consider an Alo.iIng.gSb-InSb struc- 
ture which has a potential energy depth of 100 meV. In 
Fig. [T] we show the confining potential in the longitudi- 
nal direction, Vz{z), and the eigen-functions Un{z) of the 
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single-particle Hamiltonian = -|~ Vz{z), displaced 
vertically according to their corresponding energy levels. 

En- 

The nanowhisker where the double dot structure is de- 
fined is assumed to be so thin (~ 2 nm) that only the 
lowest transverse mode is active. As such, the effective 
one-dimensional Hamiltonian of two interacting electrons 
with Rashba interaction, in the absence of a magnetic 
field, can be written aa^^ 



(1) 



where H', 







§rn^ -|- Vz{zi), m* is the conduction-band 
effective mass, zi and Z2 are the z-coordinates of the 
two electrons, and 1 and Pz,2 £^re the z-components of 
their linear momentum. Hi^r and Vint are the Rashba 
spin-orbit coupling and the electron-electron interaction 
potential, respectively. The Rashba spin-orbit coupling 
in the quasi-one-dimensional structure considered here is 
given byi^ 
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where 7n = 500 A (^xi.yi are Pauli matrices, and 



i^') ^^"^ Rashba field parameter, where the mean 



2 




80 100 

Z (nm) 



FIG. 1: The Alo.iIno.gSb-InSb double-well confining poten- 
tial in the longitudinal direction of the quasi-one-dimensional 
nanowhisker quantum dots. The single-particle eigenfunc- 
tions and energies are also shown. 



value is taken over the ground state $ of the laterally- 
confining potential Vx = Vy (T4 is assumed to be the 
same as Vy for simplicity). The electron-electron inter- 
action is given hy^ 



(3) 



Vinti\z2 ~ zi\) ^ J dxidx2dyidy2 

{xx - X2Y + (yi - 2/2)^ + (^1 - Z2Y ' 



where = (xi^yi, Zi), i = 1,2, are the electron positions, 
and e is the dielectric constant of the material (16.8 for 
InSb) . In Eq. ^ we approximate <I> by the ground state 
(spatial extent ~ 2 nm) of a harmonic oscillator potential 
defined over the cross section of the whisker. 

As a basis set for the two-electron Hilbert space we 
take 



ipi ^ ui{zi)ui{z2) |0,0), 

'P2 = -^{Ui{zi)u2{z2) +U2{zi)ui{z2))\0,0), 

<^3 = '"2(21)^2(2:2)10,0), 

ip4 ^ (mi(zi)w2(z2) - W2(2l)wi(z2)) |1, 1), 

V^S = -^{Ui{zi)u2{z2) ~ U2{zi)ui{z2))\l,-1) 

(fe = -^{Ui{zi)u2{z2) ~ U2{zi)ui{z2))\l,0). 



(4) 



The two-particle spin wave functions are the usual sin- 
glet l^) = \S ^ 0,ms = 0) and triplet states {|T+) = 
|1,1), |r°) = |1,0), |r-) = |1,-1)}. This truncated basis 



set {(pi,i = 1,6} takes into account only the two lowest 
eigenstates of ui{z) and U2iz) (see Fig. [1]). This 
truncated basis simplifies the analysis and presentation 
without leaving out essential physics. We have checked 
that the truncation introduces only small quantitative 
differences, of the order of 1 percent, in the results. 

In order to solve the eigenvalue problem of the full 
two-electron Hamiltonian given in Eq. ([1]), we expand 
the two-electron wave functions 



''ij , 



(5) 



where i — 1,..,6, and determine the coefficients by 
numerical diagonalization. It is clear that the Rashba 
coupling mixes states with different spin wave functions, 
although the mixing depends strongly in structure pa- 
rameters and applied magnetic field, as we will see in the 
next section. 
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FIG. 2: (color online) Probability density of two electrons in 
the double quantum dot with Coulomb interaction but with- 
out Rashba coupling, tpi: ground state (GS), '!/'i=2,.. .,6 excited 
states. Notice ip5 and V'e describe states with both electrons 
in the same dot, while the others can be seen more as having 
one electron in each dot. 



III. RESULTS 

We solve the eigenvalue problem of the two-electron 
Hamiltonian H given in Eq. ([T]) with the goal of un- 
derstanding the interplay between the electron-electron 
Coulomb interaction and the single-particle Rashba cou- 
pling. In order to gain some initial insight into the nature 
of the two-electron states, in Fig. [5] we show the proba- 
bility density of the six eigenstates for the two electrons 
in a double-dot structure with a 3 nm barrier without 
Rashba coupling. As expected, in the ground state the 
two electrons are essentially localized in different dots. 



3 



due to their mutual Coulomb repulsion, and on "bond- 
ing" single-particle states with a non-zero amplitude in 
the interdot barrier region In contrast, states "02 to ■04 
have electrons localized in each dot too, but with an 
"antibonding" orbital with zero amplitude in the central 
barrier. Notice also that the two higher energy states, 05 
and ipe, correspond to singlet states with two electrons 
in each dot. Furthermore, in order to fully characterize 
the two-electron system in what could be realistic exper- 
imental situations we introduce a magnetic field along 
the z-direction. The field is chosen small enough (and 
the whisker so thin) that the x-y orbital wave functions 
are not perturbed significantly by it. Thus, the field con- 
tributes only a Zeeman term to the Hamiltonian: 



J^Z — Z Oz, 



(6) 



where /i^ is the Bohr magneton, qq is the Lande fac- 
tor (go — —51 for InSb), and Sz = Si^z + <S'2,z is the 
z-component of the total spin operator. In Fig. [3] we plot 
the energy and expectation value of Sz, for all the two- 
electron eigenstates vs. the applied magnetic field with- 
out Rashba interaction, in the double dot with a 3 nm 
barrier. 
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FIG. 3: (color online) (a) Energy levels and (b) mean value 
of Sz = Si^z + 52, z vs. applied magnetic field for the two- 
electron eigenstates with Coulomb interaction and without 
Rashba coupling. Barrier width: 3 nm. \S) — singlet state, 
|r='='°) = triplet states. 



Let us point out some basic features of the results with- 
out Rashba coupling seen in Fig. [31 First notice that since 
the spatial wave functions do not depend on the magnetic 
field and only the Zeeman energy changes with B, this 
explains the linear field dependence of the energies. At 
zero magnetic field, the ground state 0i is a singlet (|0, 0) , 



Sz = 0), but around B « 2.2 T there is a level crossing 
and its spin part becomes |1, 1) {Sz = !)■ This change oc- 
curs due to the competition between the Zeeman energy, 
on one hand, and the Coulomb interaction and E2 — Ei 
on the other. At the crossing, the first excited state V'2 
goes, naturally, from Sz — itoSz — 0. 'ip2, "03 and ^4 
are degenerate at B = 0, but having Sz = 1,0,-1, re- 
spectively, their degeneracy is broken for non-zero B. E4 
and i?5 have a level crossing at i? w 6T, and, finally, 
i?5 and Eq cross at B > 7T (not shown). There are no 
further crossings at higher magnetic fields. 
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FIG. 4: (color online) (a) Energy levels and (b) mean value 
oi Sz ~ S\,z + S2,z versus applied magnetic field for the 
two-electron eigenstates including Coulomb and Rashba in- 
teractions. The strength of the Rashba coupling is given by 
i^Wx') ~ ImeV/A. Barrier width; 3 nm. I^) = singlet state, 
|T*'°) = triplet state. 



We next include the Rashba interaction with a cou- 
pling parameter (^^) = ImeV/A. The results for the 
energy levels and mean value of Sz are shown in Fig. 3] 
The main differences with Fig. [3] are as follows: 

(i) With Rashba interaction, at i3 = 0, ^02 and 04 are not 
spin eigenstates anymore, as can be seen in Fig. H^b). 
However, we still label the states as if they were pure 
spin states, as one spin state dominates the admixture 
(far from the avoided crossings discussed next). 

(ii) Two of the level crossings in Fig. [3] {Ei with E2 at 
B « 2.2 T, and £^5 with Eg at B >7T) become avoided 
crossings here, as the pair of states involved are coupled 
by the Rashba interaction. The other crossing, between 
E4 and E5 at i? « 6 T, is shifted slightly upward due to 
the effect of the Rashba interaction on each individual 
level, but it does not become avoided because the levels 
are not coupled to each other through the Rashba interac- 
tion. This lack of mixing arises from the different spatial 
symmetry of the states and the strong Coulomb interac- 
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FIG. 5: (color online) (a) Energy levels and (b) mean value 
of 5z — Si^z + S2,z versus applied magnetic field for the 
two-electron eigenstates including Coulomb and Rashba in- 
teractions. The strength of the Rashba coupling is given by 
i^Bx') ~ ImeV/A. Barrier width: 5 nm. \S) = singlet state, 
|T^'°) = triplet state, (color online) 



FIG. 6: (color online) Mean value of Sz = 5*1,2 -I- S2,z as 
a function of the Rashba parameter {^^) for the six two- 
electron eigenstates. B — 0.5 T. Barrier width: 3 nm. (a) 
Without Coulomb interaction; (b) With Coulomb interaction. 
Xso is the spin-orbit length. 



tion. The symmetry under space reversal (odd vs. even) 
prevents the mixing of a state with double dot occupancy 
(singlet -05) and a state where each dot has one electron 
("04), where each state has opposite space-reversal sym- 
metry. Notice that the width of the avoided crossings is 
determined mainly by the strength of the spin-orbit cou- 
pling, and therefore it can be adjusted with a transverse 
electric field (gate voltage). 

In Fig. [5] we present the energy levels and the mean 
value of Sz as functions of the magnetic field for a struc- 
ture with a 5 nm interdot barrier. Notice that, compared 
to the case of a 3 nm interdot barrier (Fig. [4|) , crossings 
and anti-crossings shift to lower values of the magnetic 
field. This happens because the wider barrier decreases 
the bonding-antibonding gap, the wave function overlap 
and the associated exchange integral that determines the 
singlet-triple separation at zero field. As the S-T gap is 
smaller, it is more easily overcome by the Zeeman energy. 
For an 8 nm barrier, the bonding-antibonding gap is only 
slightly over half of the value for the 3 nm barrier, and 
the singlet-triplet crossing occurs at a field of 1.2 T, for 
example. 

In Fig. [S] we present the mean value of Sz as a function 
of the Rashba parameter {^^) for all the states at a 
given magnetic field B = 0.5 T (using a barrier width of 
3 nm). In this figure, we compare two cases, (a) without 
and (b) with Coulomb interaction, in order to exhibit bet- 
ter the role of spin-orbit coupling in the spin mixing. As 
expected, without Rashba coupling, i.e. with (^^)=0, 
the spin projection of each of these states naturally (5*2) 
takes the exact values 1,0, and —1, as seen in both Figs. 



Elja) and (b). An important difference between (a) and 
(b) is that there is a symmetry around Sz = when the 
Coulomb interaction is absent. The ground state in Fig. 
[S{a) (thick cyan solid line) starts with Sz = and at a 
certain (typically large) value of the Rashba parameter 
reaches a maximum. On the other hand, in Fig.[6l^b), this 
symmetry about the zero value is lost due to the differ- 
ent mixing of two-particle orbitals ({</5i}) in higher- and 
lower-lying states produced by the Coulomb interaction. 
Figure [S] includes an axis (top) in terms of the spin-orbit 
length {Xso — fi^/'m*jR which is inversely pro- 

portional to the Rashba parameter. This length parame- 
ter helps visualize the strength of the Rashba coupling in 
comparison to the characteristic dimensions of the struc- 
ture. It is interesting to point out that spin mixing is first 
noticeable when Xso — 60 nm, the size of the two-well 
system. 

With a 5 nm barrier width, we find that the ground 
state has a mean value of Sz roughly equal to zero up to 
a Rashba constant of 2.2 meV/A in the presence of the 
Coulomb interaction, unlike the case in Fig. [S] where the 
ground state is clearly mixed for Rashba constant beyond 
1.5 meV/A. 



IV. CONCLUSIONS 

In this paper, we have investigated the effects of the 
Coulomb electron-electron interaction on the states, en- 
ergy levels, and z-projection of the spin of two electrons 
confined in a quasi-lD double-quantum-dot nanowhisker 
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in the presence of Rashba spin-orbit interaction. As a 
function of a magnetic field in the longitudinal direc- 
tion, some energy-level crossings become avoided cross- 
ings when the Rashba spin-orbit is turned on. The width 
of the avoided crossings can be controlled with a lat- 
eral gate voltage via the intensity of the Rashba param- 
eter ) . The positions of these crossings and avoided 
crossings as functions of the magnetic field can be se- 
lected by changing the width of the interdot barrier. Fi- 
nally, we found that the Coulomb interaction reduces the 
spin mixing in the ground state, as well as in the excited 



states, and displayed this reduction as a function of the 
Rashba spin-orbit length. 
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